Search for: All records

Numerical solution of partial differential equations on parallel computers using domain decomposition usually requires synchronization and communication among the processors. These operations often have a significant overhead in terms of time and energy. In this paper, we propose communication-efficient parallel algorithms for solving partial differential equations that alleviate this overhead. First, we describe an asynchronous algorithm that removes the requirement of synchronization and checks for termination in a distributed fashion while maintaining the provision to restart iterations if necessary. Then, we build on the asynchronous algorithm to propose an event-triggered communication algorithm that communicates the boundary values to neighboring processors only at certain iterations, thereby reducing the number of messages while maintaining similar accuracy of solution. We demonstrate our algorithms on a successive over-relaxation solver for the pressure Poisson equation arising from variable density incompressible multiphase flows in 3-D and show that our algorithms improve time and energy efficiency. 

The dynamics of a three-phase gas–liquid–liquid multiphase system is examined by direct numerical simulations. The system consists of a continuous liquid phase, buoyant gas bubbles, and smaller heavy drops that fall relative to the continuous liquid. The computational domain is fully periodic, and a force equal to the weight of the mixture is added to keep it in place. The governing parameters are selected so that the terminal Reynolds numbers of the bubbles and the drops are moderate; while the effect of bubble deformability is examined by changing its surface tension, the surface tension for the drops is sufficiently high so they do not deform. One bubble in a “unit cell” and eight freely interacting bubbles are examined. The dependency of the slip velocities, the velocity fluctuations, and the distribution of the dispersed phases on the volume fraction of each phase are examined. It is found that while the distribution of drops around a single bubble in a “unit cell” is uneven and depends on its deformability, the distribution of drops around freely interacting bubbles is relatively uniform for the parameters examined in this study.